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ABSTRACT:  I  wo  of  the  most  commonly  used  methods,  the  Trapezoidal  Rule  and  the 

two-step  backward  differentiation  formula,  both  have  drawbacks  when  applied  to  difficult  stiff 
problems.  The  Trapezoidal  Rule  does  not  sufficiently  damp  the  stiff  components  and  t he 
backward  differentiation  method  is  unstable  for  rapidly  varying  steps.  In  this  paper  we  show 
that  there  exists  a  one-parameter  family  of  two-steps,  second-order  one-log  methods  which  arc 
stable  for  any  test  problem  x  -  A(/)a,  Re  A {/)  <  0,  using  arbitrary  step  sequences. 
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I.  Introduction 


This  paper  is  concerned  with  the  numerical  integration  of  differential  systems, 


*),  (1.1) 

by  linear  multistep  and  one-leg  methods.  On  a  uniform  grid  (n  =  nht  n  *  0.  1 . h  >  0,  a 

linear  multistep  (MS)  method  of  step  number  ky 

k  k 

S  S  *«+  =  (1.2) 

y= o  y*o 

is  defined  by  a  set  of  constant  coefficients  |a;,  ,  y  =  0 .  k.  Throughout  this  paper  it  will 

be  assumed  that  (1.2)  is  normalized  by  the  constraint 

k 

1  (1.3) 

y-o 

k  k 

Using  the  familiar  polynomials  p(f)  =  I  and  o(?)  =  I  and  the  shift  operator 

;«n  >*o  ' 

E  defined  by  xn+1  »  £xn,  the  method  (1.2)  can  be  written  in  the  form 

p(E)xn-ho(E)f(itr  xn)  -  0,  (1.4) 


with  the  normalization 


o(E)  1  =*  1.  (1.5) 

n 

On  a  variable  grid  =  t{)  +  X  hyf  n  —  1.  2,  all  or  some  of  the  coefficients  a..  pt 

V»1  ' 

normally  depend  on  n.  For  practical  purposes,  the  M.S-mcthod  in  the  variable  step  case  is 
often  written  as 

k  k 

2  S  ^ n  +  j*  •*«+/)  K  ^  (1  *b) 

7=0  y-0 


with  the  normalization 


k 


lPj.n  = 

/  =  0 


I. 


(1.7) 


Equivalently  it  can  be  expressed  in  the  form 


(  1  K) 


2 


and  the  normalization  as 


r>„1r=o,  n.<n 

where  the  operators  />„  and  on  are  defined  by 

k 

<’„v»  =  S  a).nXn  +  /'  °nXn 

j~n 

For  theoretical  derivations  it  is  preferable  to  write  the  A/.Y-method  in  terms  of  a  step  other 
than  the  forward-most  step,  hn  +  k  =  tn  +  A  _  |.  as  is  done  in  parts  of  this  paper. 

With  every  A/S-method  one  can  associate  its  one -leg  (OL)  "twin"  or  counterpart 
which,  for  variable  steps,  is  defined  by 

ivs.-A+z/K',,-  ~  ()  n.m 

and  results  from  the  Af.Y-mcthod  (1.8)  by  permuting  f  with  the  operator  n,r  OT-methods 
were  first  introduced  in  [5|  for  theoretical  purposes  but  were  found  to  be  useful  as  integration 
methods  in  their  own  right  [  1 9). 

Our  aim  in  this  paper  is  to  look  for  methods  with  very  strong  stability  properties. 
More  specifically,  we  seek 

methods  which  produce  bounded  solutions  whenever  applied  to  the  lest  equation 

x  =  A(r ) v.  Re  \(t)  <  0  (M2) 


k 

V  p 

/-u 


(  I  10 ) 


using  any  step  sequence 

Asking  for  stability  with  arbitrary  step  .sequences  may  seem  to  be  too  strong  a  requirement. 
However,  wc  show  that  there  do  exist  methods  having  this  property  and  thus  the  above 
requirement  is  reasonable.  For  other  methods  stability  results  can  be  proved  but  only  under 
rather  complicated  assumptions  on  the  interplay  between  the  time  dependence  of  the  problem 
and  the  step  changes,  1 1 3,  20|. 

It  was  demonstrated  in  |19]  that,  if  formulas  are  implemented  in  the  multistcp  form, 
then  unfavorable  combinations  of  variable  steps  and  variable  coefficients  in  the  equations  can 
lead  to  instability.  As  a  simple  example  of  this,  consider  the  Trapezoidal  Rule: 


which  when  applied  to  (M2)  gives 


*„+l  =11+  yAn+1X(/„)lll--pirt+|X(fn  +  l)|  'v 

Let  h2rn  *  ^  aru*  *2»1+1  ®  4  and  choose  M*>  such  that  A(/2w)  =  “1  and  A(r;m+l)  =  0 
Then,  x2m  »  (-2)mjc0,  i.e.  the  Trapezoidal  Rule  is  unstable.  The  Implicit  Midpoint  Rule 
(one-leg  "twin"  of  the  Trapezoidal  Rule), 


*«+,-■*«  =  +  I  >’  |(V«  + 


on  the  other  hand  gives 


'n  +  l 


-  [1 


+  I 


+  / 


(/„  +  / 


'I  4-  I 


»l 


so  that  for  all  problems  (M2)  we  have  I  x,l+ ,  |  <  I  | 

Recall  that,  for  constant  steps,  linear  multistcp  methods  and  one-leg  methods  have 
similar  stability  properties  Jf>|  but.  as  stated  before,  with  variable  steps  this  is  no  longer  the 
case.  Thus,  although  we  sometimes  write  formulas  in  multistep  form  it  should  be  understood 
that  they  are  to  be  implemented  in  one-leg  form. 

For  constant  steps  and  constant  A’s,  our  stability  requirement  reduces  to  -stability. 
Therefore,  the  order  of  the  methods  we  arc  looking  for  cannot  exceed  two  and  it  is  natural  to 
restrict  our  attention  to  two-step  methods.  Observe  that,  for  eveiy  given  step  ratio  r .  there 
exists  a  two-paramctcr  family  of  two-step  second-order  formulas.  To  start  with,  consider 
formulas  of  this  class  whose  coefficients  depend  smoothly  on  r.  Choosing  geometric  step 
sequences  and  related  time-dependent  problems  we  show,  by  expanding  around  r  »  I,  that 

there  exists  at  most  a  one-parameter  family  of  constant-step  formulas  which  can 
smoothly  be  extended  to  variable  steps  in  such  a  way  as  to  satisfv  our  stability 
requirement. 

This  result  is  demonstrated  in  Section  5  and  was  first  given  in  |20|.  The  one-parameter  family 
of  constant-step  formulas  in  question  happens  to  be  exactly  to  set  of  those  second-order 
two-step  methods  which  are  A -contractive  in  the  max-norm,  which  in  turn  is  the  same  as  the 
set  of  methods  which  arc  /I -contractive  in  a  monotone  inner  product  norm  (it*,  (/-stable  with 
diagonal  (7) 


4 


I. cl  us  recall  this  terminology.  Wtth  any  volution  |  rn|  of  the  difference  equation  we 
associate  the  sequence  |.Y„|  where  \n  «  ( v „  ^ ^  ?  .  .  \n)  We  say  that  a  method  is 

stable  at  </  »  h\  if  it  produces  bounded  solutions  when  applied  to  the  test  equation 

x  *  A*.  AcC;  (113) 

it  is  said  to  be  contractive  at  q  =  h\  with  respect  to  a  given  norm  i|  •  II  if  !I.Yf)+j  ||  <  ||  A',,  || . 
( 1 H I.  The  set  of  stable  q's  is  called  the  stability  region  X  of  the  method  and,  similarly,  the 
contractivily  region  (  #  |j  ,s  ^ tic  set  of  <j's  at  which  the  method  is  contractive  t  he  method  is 

called  A -stable  ( A  -contractive  w  r  t,  a  given  norm)  if  the  left  half -plane  is  contained  in 
S(C  j|  ^  j| ),  and  A{)  stability  M0-conlraetive)  if  the  negative  real  axis  is  contained  in  XI C  ^  ||  ) 
For  every  A -stable  method  there  exists  an  inner  product  norm  in  which  the  method  is 
/4-contractivc,  [8|.  Note  that  when  a  one-leg  method  is  applied  to  (1  12),  then 
i*Yrt+|II  <  !|  ,Yn  #  follows  for  all  r?  and  for  any  A</)  as  long  as  |/fAf/n)|c r ^  ^ . 

Contractivily  in  the  max-norm  was  first  investigated  in  |14)  for  Adams-  and  Rungc- 
Kutta  methods  In  |l,  2|  eonlraclivilv  in  certain  polygonal  norms  was  studied  to  produce 
variable-step  stability  results  for  backward  differentiation  methods  Contractility  for  dissipa¬ 
tive  nonlinear  systems  in  inner  product  norms  ((/-stability)  was  introduced  in  |5,fi]  and  further 
analyzed  in  |7)-|ll|,  115)11*)).  Contractility  results  for  Ruuge-Kutta  methods  are  given  in 
13,4,12). 

As  mentioned  above,  the  stability  analysis  using  geometric  step  sequences  led  to 
aw* run /-step  methods  which  arc  A  contractive  in  the  max-norm  and  in  a  monotone  inner 
product  norm.  It  is  therefore  natural  to  also  look  for  variable  step  formulas  which  arc 
,4 -contractive  in  these  norms,  it  turns  out  that 

for  any  step  ratio  r ,  there  exists  a  onc-parameter  family  of  .4-contractivc  methods 

Hence  we  have  now  a)  characterized  the  set  of  constant-step  formulas  for  which  variable  step 
stability  results  might  exist,  and  b)  given  extensions  of  all  of  these  formulas  to  variable  steps 
which  actually  do  preserve  stability.  We  get  the  same  methods  using  either  the  max-norm  or 
diagonal  G-matrices.  However,  the  (7-matrix  depends  on  the  formula  parameter  and  for  every 
parameter  value  ( t.c .  any  fixed-step  formula  in  the  family)  wc  find  a  unique  variable-step 
extension  which  is  ^-contractive  in  the  same  G-norm  as  the  constant-step  formula  in  question 
Observe  that,  with  this  extension,  all  results  concerning  error  bounds  for  nonlinear  systems 
carry  over  immediately  The  max-norm,  on  the  other  hand,  docs  not  depend  on  the  method. 


and  therefore  stability  is  preserved  even  if  we  pick  an  arbitrary  ^-contractive  variable-step 
method  at  every  step.  The  max-norm  results  were  first  given  in  |20|  and  are  presented  in 
Section  3,  while  the  results  using  the  G-stability  approach  were  given  in  [  1 0|.  and  are 
presented  in  Section  2. 

In  Section  4  we  analyze  ,40-contractivity  for  variable  steps  and  give  variable-step 
extensions  for  all  constant-step  formulas  which  are  /t0-contractive  in  the  max-norm.  We  also 
briefly  discuss  the  selection  of  particular  methods  from  among  the  one-parameter  family  of 
/4-contractive  methods.  Finally,  we  rewrite  all  methods  in  a  practical  parametri/ation  and  list 
their  error  constants. 


2.  /4-contractivity  in  the  G’-norm 


Consider  the  general  variable-step  k- step  formula 

Srt/  rn  +  7  =  h"  £  P, 

,r-A\  /-0 


with  the  normalization 


* 


S/<,= ' 


(2.1) 


(2.2) 


For  k  >  2*  the  coefficients  ar  fit  depend  on  the  step  ratios 

r»*j  =  h„n/h"+J~  f  >  "  2-  •  *• 


(2.3) 


where  /i/f  =  tn—tn_  ],  although  for  simplicity  in  writing  they  are  not  subscripted  here,  and  hn  is 

some  homogeneous  function  of  the  first  degree  of  the  hn±r  j  1 . k ,  The  formula  (2.1) 

can  be  written  using  the  difference  operators 


k  a 

px, .  :=  S  **/  •  °v«  :=  S  v»+/ 

/  =  0  /  -  u 


(2.4) 


where  p  =  pM,  a  =  In  terms  of  these  operators,  (2.1)  becomes 


PXn  =  h'OX„.  (2.5) 

Two  possible  implementations  of  (2.?),  as  applied  to 

x  =  /(r,  x).  (2.f>) 


are  the  familiar  linear  multistcp  method 

(”fn  =  *„"/<'„•  *„> 

on  the  one  hand  and  the  one -leg  method  15,6] 


PX„  =  hJ(otn,  «,) 


(2.7) 


(2.8) 


on  the  other  hand.  Necessary  and  sufficient  conditions  for  the  onc-lcg  method  to  be  of  order 


of  accuracy  p  were  given  in  [l()|  in  terms  of  the  moments 

k  k 

MAp)  :*  £  Vr  :®  X  /Vr  *'  ®  °'  1 . 

7-0  7-0 

where  /  »  /*  +  j  =  0 . A:,  and  /*  is  some  reference  point.  I  he  method  (2 

order  p  iff,  in  addition  to  the  familiar  linear  constraints 

A/r(p)  «  *Mr_  j(^).  »'  =  0,  I.  .  p 

for  />lh  order  accuracy  of  the  mullistcp  method  (2.7),  the  nonlinear  conditions 

-  l m r  1  *  *  i . r 

are  satisfied.  For  p  =  1 ,  (2.10)  for  r  =  0  and  I  represents  the  consistency  conditions 

k 

,,m=  = 

7-0 


S-v,  =  "(<)=  '• 

»-<»  /»0 

the  second  of  which  coincides  with  (2.1  U  for  r  =•  I  in  this*  ease.  Also,  for  r  =  2.  (2 
identity  and  thus  the  condition 

M=0  7«t> 

added  to  (2.12)  guarantees  second-order  accuracy  for  both  (2.7)  and  (2.X) 

Consider  now  the  general  two-step  formula,  i  c.  (2.1)  for  k  ==  2  with  / 
T,  *  r,  n  =  0,  and  with  r0  *  r0  t2  «  t2  „  defined  by 


hn*\  =^.|-^  =  -To/S- 

+  =  TA 


(2.')) 


.X)  is  of 


(2  10) 


(2.11) 


(2  12) 


1  1)  is  an 


(2  1  31 


M  ~  ;M+  1* 


(2  14) 


In  this  cast*  l he  consistency  conditions  (2.12)  become 


l>(  I  )  —  d 1 1  4-  (t  j  +*  n  =  0. 


(2  15) 


In  addition,  the  constraint 


'(!);=  V  <t(U 


is  being  imposed  for  the  purpose  of  pinning  down  htr  and  ::  (see  (2  24)  hereafter). 
Relation  (2.1b)  is  reminiscent  of  the  first-order  consistency  condition  (2.12b)  but  is  equivalent 
to  that  condition  only  in  the  case  of  equal  steps,  re  when  hr  4j  =  =  ht.  =  h  and 

-To  =  T:  =  1 

From  (2  15)  ami  (2.16)  it  follows  that  n(>)  >>(0  =  «(.*>  m.-)'  lor  : 

(f  +  ))/(;-!)-►»  (it*  as  >  1  );  here  Mr  I  (.•  lYvit  <:)).  '(.->  :=  (.*  -  I  )'<*(((:)).  .irul 

$(r)  =  (r  4  l)/(r—  I).  Ilcncc  there  exist  i|iiantiiics  <r.  />,  ami  <  such  that 


2 .’  ( i  ( *  >  >  /’(( (r))  -  r  4 


rt:  =  1  (.•  4  I)  .  /»,  =  1  |  1  4  h  4  (.J  4  •  )|. 

2  4 


n  i  =  -<  * 


( 2  1 S ) 


nn-  i 


P„«  -!  1 1+ *-<«/  +  <■>! 

4 


By  (2.12),  (2.13)  and  (2.  IX)  the  first-  and  second-order  accuracy  constraints  can  be  rewritten 


(r  +  l)r^  +  ( e—  1)t„  -  2, 


(2  1<>) 


(c  4-  l)r;  +  (c—1)tJJ  =  (I  +  b)(r2  +  T0)  f  (<;  +  «  )(:,-:,,) 


(2  20) 


9 


The  variability  of  the  step  si/e  can  he  represented  by  the  parameter 

o  :=  <•( t0  4  t:)/2  (2.21  ) 

which  vanishes  for  equal  steps.  From  (2.19)  and  (2.21)  we  get 

T  2  ~  T0  =  2(l-o).  (2.22) 

The  second-order  accuracy  condition  (2.20),  which  is  equivalent  to 
c[ ( t2  +  T0)fc  4  (t2-t0)a1  +  2(t2  4  t0)(t2-t0)  =  2(1  4-  b)(r2  +  r{))  4  2 (a  4-  (*)(t2-t0), 

can  thus  be  rewritten  as 

(*[  (4o *  / )  +  4 (  1  —  cr )  ^  |  4-  X (  1  —  a  ) a  / C  =  4(1  4-  h)n/c  4-  4(0  4-  C ) (  1  —  a  ) . 


or 


(a  —  rt<’)(l—  f*)  =  bn  4-  or(l  —  or). 

For  a  =  0  (equal  steps)  and  c^O,  this  last  condition  reduces  to  a  =  0,  and  for  n(l-n)/0  it 
can  be  written  in  the  form 


or 

a 


+ 


b 

1  —  n 


1  -c\ 


(2.23) 


From  (2.14),  (2.21)  and  (2.22)  it  follows  that 


=  1 -a  4  j. 


-T() 


l  —  o 


r 


(2.24) 


/i 


»» 


=  \.  +  2/t2  * 


The  formal  algebraic  condition  for  A-stability. 

Re  p(? )/o(f)  >  0  for  |f  |  >  1. 


is  equivalent  to 


/(r  r(:)/,i(;)  >  0  for  /lr;  >  0, 


i.e.  that  r{z)/s{z)  be  a  positive  function  |9|;  here  r(*)  =  2(:  4  c)  and  a(c)  - 
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jU  +  c)  +  az  +  b.  From  (2  17)  it  follows  that  this  condition  is  satisfied  iff 

a  >0.  b  >  0.  e  >  O  (2.25) 

An  interpretation  of  the  condition  (2.25)  as  an  actual  (not  just  formal)  stability  result  is  given 
in  Section  5  hereafter.  It  is  useful  in  the  following  to  carry  out  the  investigation  of 
/f-contractivity  subject  to  this  constraint. 

To  analyze  /f-contractivity  of  (2.1)  in  the  (/-norm,  recall  |7)  that  the  determination  of 
a  (7-matrix  is  related  to  decompositions  of  the  form 

Re  jr(;|i(2)  =  (  |  m,r  +  wi0  |  ‘  +  n„).r  +  |C|(;)|‘  +  t  vt,<- >  I  \  (2.20 

where  x  —  Re  z%  nt  |  >  0,  nu  >  0,  and  ^ ,  are  linear  functions.  In  fact. 


where 

mo  +  »i 
"*|Wln 

Note  that  G  is  positive  definite.  In  homogenous  coordinates  P  corresponds  to  the  Mocbius 
transform  z  as  (f  +  !)/(£*-  I)  Put  for  brevity, 

*n  *i:l 

*12  *22  J 


7  A 

G  =  PlGP , 


(2.27) 


?]  •!  :] 


We  may  write  (2.26)  in  the  form 

Re  ±r(zw7)  =  (1.5.  52)(1g  +  //)(!,  .rV  (2.2S) 

2  2 

where  H  is  a  positive  semi-definitc  matrix,  the  elements  of  which  arc  determined  by  the 
polynomials  ^ ^0. 


*11  *12-1  p *  1 1  h  12  °-i 

*11  2*I2  *22  .  W  -  I  A, 2  *22  <>  I  >  0- 

*12  *22  0  0  0  0 


One  decomposition  of  the  form  (2.20)  is  easily  found,  namely 


1 1 


Re[^r(z)x(2)\=  Re\(z  +  c)\}(i  +  < )  +  az  +  b\\  -  |  | ;  +  < |  *  +  (<i<  +  b)  |*  +  a  \  :  |  *  +  he 
In  this  .special  case  we  obtain. 


We  are  interested  in  all  possible  decompositions  of  the  form  (2.26)  or,  equivalently.  ( 2  2 S ) .  It 
follows  from  (2.2S)  and  the  form  of  G  and  //  that  /i||,  gn  +  2/ii:.ir,7  ef<*.  are  the  same  for 
all  such  decompositions.  Hence. 

/i ,  |  =  be,  g 1 2  ~  c,  =  I,  ss  </.  2/i , 2  =  e"  +  </c  4*  /*-,e ,, 

A 

In  addition  to  this,  we  have  the  positivity  conditions  for  (/  and  //  Put  for  brevity  =  g 
Then 


<7 


e  > 


( 2. 2(i ) 


and 


d/v 


>  ^i: 


i(r 


+  II,  +  /•-cl* 


(2.30) 


A 

The  matrix  (V  associated  wilh  (/  of  (2.21))  is 


G  - 


(1  +  ,e-2<*  l-.e  \ 
I  -g  \lg  +  2c) 


(2  31) 


it  is  subject  to  the  constraints  g  >  c*  and  (2  JO) 

Global  ,4-contractivity  (and  thus  .4-stability )  results  are  obtained,  for  all  solutions  of 
monotone  nonlinear  systems  |6]  generated  by  vat  table-step  nnc-lep  methods  of  type  (2.K),  if 
for  all  n  >  0  the  formula  (2.1)  is  ^-contractive  in  a  (/-norm  with  a  G  which  is  inJcpctulenf  of  n 
(i.e.  of  the  step  ratio).  For  G  of  (2.31)  to  be  independent  of  w,  g  and  c  must  be  constant. 
The  latter  means  that  the  difference  operator  p  defined  by  (2.4)  is  constant  with  respect  to  n 
From  this  point  on  the  analysis  will  be  restricted  to  second-order  formulas,  f  c. 
formulas  for  which  (2.23)  holds.  Since  in  this  ease  a  0  if  n  0.  (2.3f»j  implies  that 
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0)  =  0  and  that 


b(a )  +  a{a)r  —  g  -f  /•*  «  0(  1  ain )  |  ‘  J.  n  -►  0. 


Hence  />( )  is  continuous  at  nr  —  0.  From  (2.23)  it  follows  that 


MO)  —  1  -*r‘  —c  lint  (a / 


i.e,  lint  {a/a)  exists.  Since  a  >  0  for  |a|<<  1  because  of  the  .^-stability  requirements,  it 
o— ►o 

follows  that  lim  (a/a)  =  0  (since  the  right  limit  is  >  ft  and  the  left  limit  is  <  0)  Hence 
2 

b{ 0)  =  l— c  ,  g  =  1,  and  the  only  possible  choice  for  G  is 


G  =  2 


C; • 


<  V  <  I, 


(2.32) 


a  diagonal  matrix 


If  in  (2.30)  one  lets  g  =  I  and  uses  the  second-order  accuracy  condition  (2.23)  to 


express  ac  in  terms  of  h.  a,  and  c  then,  in  this  cn^c,  (2.30)  becomes 


—  |  (  I  —  a )( I  —  rM )  —  —  |  >  ft 

I  —  a 


(2.33) 


Therefore,  the  onlv  second-order  methods  which  are  A -contractive  in  a  fixed  G-norm  for 


arbitrary  step  ratios  are  defined  by 


a(« )  =(l— c")a2/r, 

hi  a)  =  ( 1  —  c  " )  (  I  —  «)  . 

r(n)  a  r  a=  const.  0  <  c  <  1. 


(2.34) 


Another  way  of  arriving  at  the  one-parameter  class  of  methods  (2.34)  is  to  observe  [10|  that 
for  any  fixed  c  the  equation  (2.23)  defines  a  onc-paramctcr  family  of  straight  lines  in  the 
plane  with  cartesian  coordinates  (ac,  b).  The  envelope  of  this  family  is  given  by  (2.34).  As 
mentioned  above,  in  order  to  have  a  G-norm  independent  of  w,  the  difference  operator  p  must 
be  fixed. 


*  \ 


3.  /l-contractivity  in  the  ^x-n«rm 


In  the  preceding  section  it  was  shown  that  the  class  of  all  second-order  two-step 
formulas  with  variable  steps  which  are  A  -contractive  in  the  G-norm  is  defined  by  (2.34).  Here 
it  will  be  proved  that  this  same  class  also  represents  all  variable-step  p  =  k  =  2  formulas 
which  are  /t-contraclive  in  the  ^-norm  However,  as  will  be  seen,  the  requirement  that  the 

operator  p  be  constant  with  respect  to  n  which  had  to  be  imposed  in  analyzing  /f-contractivitv 

% 

in  the  G-norm  may  be  dropped  in  studying  the  same  property  in  the  — norm. 

Let  h~  :=  /t„  +  1  =  h+  :=  A„  +  ,  =  +  ami  h  :=  i(/j+  +  ID.  Then 

one  has  h +  =  /i(l  +  f)  and  h~  =  A(  1  —  r ).  where  f  =  (/i  +  —h~  )/2/i,  -1  <  r  <  +  I:  /  e 
c  =  0  for  equal  steps  =  h*).  For  bounded  step  ratios,  one  has  the  following  expressions 
for  the  first  and  second  central  divided  differences  of  any  sufficiently  smooth  function  ,v(r): 

- ^ — |  —  ( l  +  r)2v(/„ )  +  4rx(tn  +  , )  +  ( I  -f  )‘x(»n  +  ,)|  =  )  +  (>lA‘). 

2(1  -t)h 

(3  I) 

- L — -|(  l  +  r).x(r„1-2<(f,(  +  l>  +  ( 1  —  r ) v(fB  + 1)|  =  \V,,4 ,  )  +  0(/i). 

(l-f  2)lr 


If  unsubscripted  quantities  are  assigned  to  /M4_r  then  one  can  write  a  notmali/ed  ansatz  for  a 
two-step  formula  in  divided  difference  form: 


.fix  ,  lt2fi* 

«<><  +  u\h  77  +  uzh  - 

Of  X, 


-  -  /if  x  +  v ,  h  -V- 
2  V  1  fit 


In  order  for  (3.2)  to  be  consistent  (of  order  of  accuracy  p  =•  1 )  it  must  be  satisfied  exactly  for 
*(/)=!  and  *(/)  =  t.  This  is  the  case  iff  w0  =  0  and  //,  =  I.  The  class  of  all  consistent 
formulas  can  thus  be  written  in  terms  of  the  three  parameters  u  =  u2,  v  =  i  , .  and  w  =  v2  as 
follows: 


h--  +  ltlr  - h(  X  +  vh  +  H -h2-'*  )  =  0.  (3.3) 

f>t  k,z  ' 


The  p  »  2  accuracy  condition,  gotten  from  substituting  v(/)  =  /2  into  (3  3).  is  simply 


i/  =  v.  (3  4) 


If  from  (3.1)  one  substitutes  for  the  divided  differences  and  groups  terms  one  finds  the 


following  expressions  for  the  formula  coefficients 


«,)  *  Y -A. 
<»1  =  ~2>* 
<>2  =  Y  +  <T 


/t((  -  , 

/*,  *  1-20. 

/f  ,  =*  0  4-  v'  . 


where 


Y  •* 


-2#  -f  2u 
2(  I  — r2  ) 


I  4  r  ‘  —  2  nt 
2(  I  ~  f : ) 


—  2r  i  4  2m 
2(l-r:)  ’ 


(  I  -f  f  ‘  )v—  2r  ve 
2(  l — r 2 ) 


Two  identities  are  implied  by  these  definitions 


I  -  2#  y  -2<5. 


i  0  +  4  —  - 
2 

Note  that  the  normalization  is  again  (2  2).  i.e.  1  /(  =  I 

/»«  ' 

The  class  of  formulas  defined  by  (3.5)  and  (3  6)  N  the  same  as  that  defined 
of  the  parameters  a ,  />,  c  by  (2  IX)  The  transformation  ( «/ .  v.  h  )-+(<;./’,  <-)  and 
related  quantities  can  be  found  by  requiring  that  n  f  h  =ou  x  M  and  h  I 

h  =  /in.  One  obtains  the  relations 

a  =  4  ^  —  r  =  4  v'  —  ( Y/<^ ) . 

/)  =  4</>  —  1 , 

c  —  y  t  fi, 
a  =  2ry  =  I  —  2^ , 

h  =  A/2.S. 

r  j  =  2(1  +  r)<S. 

-r<|=s  2(l-f)A 

In  the  present  paramclrization,  r(z)  =  4(fo  4*  y)  and  \{:)  =  ::  4  4^:  4  Mis¬ 

applying  the  result  of  Fxamplc  2.  p.  1 3  of  |9|  to  iU)/rO)  one  finds  the  formal  A 


(3  5) 


( 3  f. ) 


(3  7) 


in  terms 
that  of 
fi.  where 


(3  X) 


I  )  IU 
stability 
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conditions 


A  >  0. 

<t>  >  1 .  (.■>')) 

4 

0  <  y  <  4  fill 

It  was  proved  in  jlH|  that  a  multistep  formula  is  .4,, -contractive  in  the  max  norm 
(^-norm)  iff 

*k  >  0.  «1;  <  0.  /  *0.  ..  k~  I. 

(3  1(0 

k  i 

ih  -  V  l  il  |  >  i) 

fall 


A  formula  is  ^-contractive  in  the  max  norm  iff  both  (3  10)  ami 

k  •  i 

W/Ci|)  <  I.  0  <  r  )<+/-.  (.Vi  1 ) 

,*(l 

arc  satisfied,  where 

K;tn)  »  [(n‘  V  /V  q  +  /^rj)J  .  /  *  0,  ...  A  -  !.  (  V  1  ^  > 

and  where  the  square  roots  are  positive  by  definition  I  or  A  =  2  it  is  le.isible  to  reduce  t  >  II) 
to  a  condition  just  involving  the  formula  coefficients  In  fact.  Rn  +  R  ,  <  I  <:-•> 
2A?0/( ,  <  4^  4 R^R]  <  1 1  ~*KC*;— At  J the  last  inequality  being  rational  1  hu 

latter  is  in  turn  equivalent  to  a  condition  of  the  form  p{)  4*  p,T|  4  poi*  >  0.  17  >  0.  m  which 
p()  *  0  by  consistent')  and  p,  >  0  by  (3.10)  therefore,  a  two-step  formula  u  ,f -contractu  c 
in  the  max  norm  iff  both  (3.10)  and 

p,  tn?  -  nj  4  -  /if  +  /if, »  -  2l. 1 4  >0  (I  M) 

hold  After  substitution  from  (3.5).  a  calculation  gives  the  condition 

«  (4'->‘)(4^-  I )  —  4 ( y ^  >  (l  <3  14) 

which  is  equivalent  to  (3.13).  Using  the  notations  A  : -  -*/(!-#").  R  :=  (l-r‘)  \  and 
C  :»  (I  +  r2)/2(l-f")  and  the  second-order  accuracy  condition  (3  4 ) .  one  has  y  =  A  4  /n. 


6  *  C‘  +  4>  *  Av  +  B h ,  and  ^  *  (V  • f  -Im  Anothei  calculation  vie  Ids  - 

~  —  Ax  —  Bv2,  4<>— I  *  4 Av  +  4fln  -  I,  4(yif  — »V.)~  =  | /? ( ***  —  »~)|  .  and  thus  (3.14)  finally 

4 

takes  the  form 

1,,,  =  -  t  |(  |  -4.I*)-2«(h  +  »:||‘  >  0 

This  last  condition  is  satisfied  iff  *  =  ( 2B )  1  -{2A  ' B)v~ v\  t  e  iff 

H-  -»  +2(»-t*’=  +'  (.VI5.I) 

2  ^ 

As  will  he  shown  in  the  next  section,  <3  10)  implies  that  ol  the  parabola  (3  15a)  only  the  are 
defined  by 


f  <  v  < 


I  +  f 
2 


1  <i  <  I . 


(3  15b) 


corresponds  to  ^-contractive  methods  It  can  be  verified,  hv  using  the  transformation  (3.8), 
that  the  class  of  second-order  formulas  defined  by  (2.34)  is  the  same  as  that  defined  by  (3  3). 
(3  4)  and  (3  15). 

Remark  The  use  of  the  f  r-norm  is  natural  since  il  is  the  only  norm  in  which  all  of 


the  formulas  defined  by  the  interval  (3.15b),  incliiifinr  its  endpoints  the  Jrape/oidal  Rule 
{v  =  (I  +  t)/2)  and  the  two-step  Trapezoidal  Rule  (v  =  r),  are  simultaneously  A -contractive 
( 1 20),  Proposition  2  4). 
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4.  40-con*ractivify  and  summary  of  methods 

It  was  proved  in  [18]  that,  in  the  ^-norm,  the  formula  (2,1)  is  a)  contractive  at  q  - 

0  iff 


<*j  <  0,  j  =  0,  ..  *  k—  I ,  >  0, 


(4  1) 


b)  contractive  at  q  =  *  iff 

y  =  <'• 

./*  0 


(42) 


and  c)  /40-contractivc  iff  both  (4.1)  and  (4.2)  hold  (as  slated  earlier)  For  the  class  of 
formulas  with  p  =  k  *  2  defined  by  (3.4)-(3.6).  n0  <  0  and  o,  <  0  yield 

f  <  »■  <  .  (4.3) 


The  condition  >  0,  which  is  equivalent  to  v  >-(  l-r)/2.  is  implied  by  (4.3)  In  discussing 
(4.2)  one  distinguishes  four  cases  depending  on  the  signs  of  (<0  and  /(,.  first  note  that  0 
iff 


> 

H* 

< 


Llr 
2  ‘ 


(44) 


and  /?,*  0  iff 
< 


Cave  l :  For  /?n  >  0  and  /( ,  >0  one  has  y  =  ~  ^  +  v'  )-*  I  >  0  iff 


H  >  — 


1  +  r 


l-_r 

"  2 


Case  2:  For  ft{)  <  0  and  /?,  >  0,  y  *  /?2  +  /*n-/*i  *  4<f>-  *  £  0  iff 

M»  >  —  ”  f  +  f  l\ 

4 


(4  b) 


(4  7  | 


l-r* 


4 


+  n’  < 


I  4  t 
~  2 


But  (4.7)  and  (4.4)  require  that 


IS 

which  implies  i  >  ( I  +  t)'2  I  his  latter  condition  is  incompatible  with  (4  1),  so  that  Case  2 
does  contribute  any  feasible  points  to  the  region  of  .•f0-contractivity  in  the  parameter  plane 
Case  .1.  I'or  /I,  <0  and  /?0  >  0.  y  =  +  2( V  -<M  +  1  >  0.  /  e. 


w  < 


I  +  r 


v  + 


I 


( 4  K ) 


l' use  4.  fair  /fj  <0  and  /f„  <  0,  an  argument  similar  to  that  of  Case  2  shows  that  there  are 
no  , -^-contractive  methods  of  this  type  I‘he  region  of  .f()-conlractivity  in  the  plane  of  the 
parameters  (n*t  v)  is  the  shaded  triangle  shown  in  f  igure  1  for  increasing  steps  {r  =  ■— )  and  in 
f  igure  2  for  decreasing  steps  (f  =  )  The  line  -0  divides  this  triangle  into  two 

subtrmngles.  the  left  hand  one  represents  the  contribution  from  Case  I,  the  right  hand  one 
that  from  Case  3  1  he  upper  left  anti  lower  vertices  of  the  4,rconlractiviiy  triangles  represent 

(identically  in  r).  ihe  Trapezoidal  Rule  (TR).  defined  by  —  /f0  =  0  and  the  two-step 
Trapezoidal  Rule  (21R).  having  « ,  =  /t(  =0.  The  set  of  all  second-order  A -contractive 
methods  (both  in  the  G-norni  and  in  the  /  —  norm),  defined  by  (2  34)  respectively  by  the 


cquiv;dcnt  relations  (3  15),  is  reprrsentcil  by  ares  of  parabolas  in  I  igures  I  and  2.  joining  TR 

and  21  R 

One  may  ask  whether  it  is  possible  to  select  a  particular  formula  from  among  the 
,4(l-contraclive  or  the  d-contractive  ones  whose  truncation  error  is  in  some  sense  minimal?  Of 
course,  there  is  no  unique  way  of  making  such  a  choice  If  the  problem  is  smooth  and  no 
dissipation  is  required,  the  TR  (resp.  its  one-leg  countcr[>art,  the  Implicit  Midpoint  Rule)  or  a 
method  "dose”  to  it  may  be  the  tight  choice,  while  for  rough  problems  one  needs  a  considera¬ 
ble  amount  of  dissipation  which  one  can  provide  by  asking  for  strong  contractivity  at  q  ~ 

In  |ll)|  a  compromise  choice  of  an  40-contractive  formula  was  made  b\  minimizing,  over  all 
At) -contractive  formulas  with  p  =  A;  =2.  a  bound  for  the  global  truncation  error  produced  if 
an  ,4(0)-eontractive  formula  is  implemented  as  a  one-leg  method  and  applied  to  the  test 
problem  i  =■  M / )\,  \(t)  <—  u,  a  >  0.  Ihe  calculation  was  restricted  to  fixed  integration 
steps  It  resulted  in  finding  an  4(, -contractive  (indeed.  -Hal-contractive)  formula  of  Adams 


type. 


r„.  i-*<  l « 


n  +  l 


O  -  0. 


(43)) 


which  in  [ll)|  was  referred  to  as  CA2  A  unique  generalization  of  (4  <))  to  arbitrary  variable 
steps  was  then  defined  by  requiring  that  the  variable-step  version  represent  a  seeond  order 
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accurate  Hermite  interpolation  formula  of  the  same  form  as  (4  9),  /.<■  with  n(l  =  /*,  =  (),  and 

it  was  proved  that  the  variable-step  version  defined  in  this  manner,  when  implemented  as  a 

one-leg  method,  remained  /f0-contractivc  (in  the  max  norm)  for  any  A(/ )  <  0  and  for  any  step 

sequence  whatsoever.  For  a  given  k  the  objective  function,  whose  minimization  over  the 

two-parameter  set  of  all  .40-contractive  formulas  led  to  (4.9).  was  |c,  |  y;  here  c a  is  the  error 

constant  and  y  =  0»  —  2  \fl  \  is  the  contractivity  constant  considered  in  (3.10)  and  (4.2). 
7«<>  J 

both  of  which  are  functions  of  the  parameters.  It  can  be  shown  -  although  it  was  not  stated  in 
1 19]  -  that  formally  minimizing  It^l/y  over  all  /fw-eontraetive  formulas  in  the  variable  step 
case  defines  the  same  generalization  of  (4.9)  as  does  the  Hermite  interpolation  requirement 
mentioned  before.  The  variable  step  version  of  (4.9)  is  defined  by  (4.19)  hereafter  in  the 
form  given  in  ]19)  which  is  practical  for  numerical  computation.  It  is  represented  in  Figures  I 
and  2  by  the  point  labeled  OPA{)  (for  "optimal  ”  d^-conlractivc  method),  the  intersection  of 
the  lines  «0  »  p{  =  0. 

The  transformation 


M.(f)=  U  +f)l  ( I  +  r  *  )h’(0)  —  f  "vtO)  1 

(l-r)  +  2r|2>v(0)~v(0)| 

(4  10) 

_  (1  -f  )|2in(())  +  (l -2r)i(())l 

( 1  —  r  )  +  2t  1 2  w(())  — 1*(0)  | 


maps  the  triangle  representing  the  /f0-conlraetivc  formulas  for  f  =  0  one-to-one  onto  the 
corresponding  triangle  for  any  f*0.  The  lines  n()  =  =  /*,  =  0  are  invariants  of  this 

transformation  and  thus  the  fixed-  and  variable-step  versions  of  OPA{]  and  of  the  backward 
differentiation  formula  (Bl)F)  correspond  to  one  another 

Wc  derived  a  particular  ,4-conlractivc  />  «  k  =  2  formula  for  arbitrary  step  ratio  which 
minimizes  the  objective  function  |e,|/y  mentioned  in  the  preceding  paragraph.  I  he  error 
constant  for  the  p  =  k  =  2  formula  associated  with  (v,  h  )  turns  out  to  be 

c,  =  [ ( I  —  r ‘ )  +  4rr-ftn-];  (4  11) 

on  the  class  (3.15)  one  finds  that  |r,|  =  ^-[(1  +  *")  +  4m-3v"J  A  calculation  shows 
that  |  Cy(\)  |  /y(v)  takes  its  global  minimum  at 

V,vl  m  (1  +  r):(3  +  f)"1  (4,12n) 

which  is  an  interior  point  of  (he  interval  defined  by  (3.15b)  for  all  r.— 1  <  r  <  1.  I  he  value 
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of  H\  corresponding  to  vopf  via  (3.15a),  is 

*',ri  «  —  ( I  +  f  )(7  +  +  5#  “  +  <  ')  '(3  +  f  )\  (4. 1  2M 

For  equal  steps  (r  =  0).  the  particular  A  -contractive  formula  defined  by  (4  12)  is 


Ir-4, 


*  +  ! 


+  +  i'”"  +  ip-*?)  =°  (4 13) 


It  is  represented  in  figures  1  and  2  hy  the  points  labeled  OP  A  for  <  "optimal”  /(-contractive 
method) 

For  practical  use  the  family  of  variable-step  p  —  k  —  2  methods  studied  in  Sections  2 
and  3  is  rewritten  hereafter,  with  n  —  2  replacing  in  terms  of  the  forward  step  hr  =  1 

(rather  than  the  average  step  h  =  used  in  Section  3)  and  in  terms  of  the  step 

ratio  r  m  rn~  hn/hn^  (rather  than  r  =  (hn  —  hn]  )/2h\.  f  he  relationship  between  r  and  r  is 


I  +  f 

T-r 


r-\ 

r  +T 


(4  14) 


The  class  of  formulas  defined  by  (3  5)  and  (3.(0  (or  (2.  IS))  is  equivalent  to 


2  +  +  P\xn-\  f  P:x^  = 


(4.15) 


where  the  coefficients 

«o=~i£;<‘-V  •  ft,- 

«,  =  —(!  — r(  I  —/>,)]  .  /i,  =(|-1/)||-^(1  +  />„-/»,).  (4.l(i) 

„  =1-  r_(|_fr)  ,  f]  m  1/,  + - r - (1  +  /*„-/», ). 

2  I  +r  1  2  2  1  2(l+r)  0  1 


depend  on  n  via  the  step  ratio  r  *  rn  and  possibly  via  I  he  parameters  />0,  />(  which  may  also 
vary  with  n.  The  parameters  v  and  h  of  Section  3  relate  lo  h{),  h ,  as  follows: 


l+r  . 

v=  ~Y~h\- 


H’  *  (  )?(  I  +/>0) 


(4  17) 
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The  use  of  b0  and  b{  is  motivated  by  the  fact  that  for  r  *  1  (equal  steps)  one  has 
s(z)  =  o(z(f))  =  b{)  +  bxz  +  2 and  the  /(-stability  set  in  the  parameter  plane  (/>u.  />,)  is  the 
first  quadrant.  The  error  constant  for  the  general  formula  (4.16)  is 


Cj  =  -p-— 1(2  — 3r)  — 3rA0— (2  — 2r)fc,  |. 


(4. IS) 


The  variable-step  version  of  the  optimal  /f0-contractive  formula  (4.10)  labeled  OPA()  for 
which  b{)  =  1/r  and  bx  =  1 ,  is 


xm  —  xH 


i-a„[ 


2  +  r 
2  +  2  r 


2  +  2  r 


*-»]  - 


(4  19) 


and  its  error  constant  is 


<\v 


Hr 
1 2r 


The  constraint  (3.15)  (respectively  (2.34))  for  /l-contractivity  iranslales  into 


h{)  = 


r 


r—  1 
r 


<  *1 


<  I 


(4.20) 


(4.21) 


The  a-cocfficicnts  of  the  /(-contractive  formulas  are  those  of  (4.16)  and  the  /^-coefficients  are 


P  0  = 

p,  = 

— : —  + 

l+r 

l~rb{ 

27urT  1 

2(  l  +r)  '' 

(4.22) 

r  1 

07  = 

2  1 

1  . 

2 

2r-l  h 

r  h2 

l+r 

2(l+r)  1 

2(  1  +  r) 

The  error  constant  for  the  /(-contractive  methods  is 


Cy  *=  -i-f-4  +  4( I  — r)/> |  4  3r/)J] 


(4.23) 


A  variable-step  extension  of  the  formula  OVA  given  by  (4.13)  is  defined  naturally  by  the 
relation  (4.12a),  i.e.  by  requiring  that  the  variable-step  version  have  the  same  minimality 
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property  as  the  uniform-step  version.  Condition  (-1  12a)  translates  into 


/>,  =  2r  ( 2r  +  |). 


the  formula  coefficients  are 


(  l+r)(l+2r)  * 
l+r 


1+2  r+2r 

2  ~  ( 1  +  r )  ( 1  +  2r )  ’ 


/*o  =  /<. 


r(  I  +r) 

( i  +:d: 


Pi 


I  +2r+2r: 
(1+2  r)2 


and  the  error  constant  is 

__  _  (1  +r)(J  +  ^+C) 
3r(j+2r): 


(4.24) 


(4.25) 


(4  2(>) 


Another  extension  of  (4.13)  is  given  by  the  requirement  that  the  variable-step  version 
be  A  -contractive  in  the  same  (7-norm  as  (4.13).  Using  the  notations  of  Section  3,  this  requires 
that  v  =  v(r)  =  (l  +  3r  +  rA)/(3  +  2r),  resp 

bAr)  =  (4.27) 

r<5r+|) 

Recall  that,  when  formulas  are  implemented  in  the  one-leg  form,  the  local  truncation 
error  consists  of  two  third-order  terms.  The  first  of  these  has  the  coefficient  c,  and  comes 
from  the  linear  constraints  (2.10);  the  second  term  is  due  to  the  nonlinear  constraints  (2.1  1 ) 
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5.  On  variable-step  stability 

As  stated  in  the  Introduction,  wc  seek  extensions  of  uniform-step  p  =  k  ~  2  formulas 
to  arbitrary  non-uniform  steps  whose  coefficients  a  lif  depend  smoothly  on  f  near  t  =  0 
(uniform  steps)  and  which  generate  bounded  discrete  solutions  when  applied  to  the  variable 
coefficient  test  equation  (1.12).  A  necessary  condition  for  a  method  to  have  this  property  is 
that  it  be  stable  for  any  geometric  step  sequence 

hn  =r  (1  4*  r)”  1 ,  f  constant.  (5.1) 

for  which  the  step  ratio  is  constant: 

r  =  1  +  r .  f  =  r/(2  +  r).  (5.2) 

and  thus  the  operators  p  and  o  are  constant.  Furthermore  let  A ( / )  be  such  that 

\(otn)  »  (1  +  r)1  nq,  q  constant.  Re  q  <  0.  (5.3) 

n 

where,  for  t{)  =  0,  tn  =  X  hv  =  [(1  +  e)  -!)/>.  Then  q  =  =  constant  and  the 

Vss  1 

one-leg  method  generates  bounded  solutions  iff  p( f,  f  )/*(£.  r)  satisfies 

I  f  I  >  1  =>  /te|p(f.  f)/o(?.  r)|  >  0.  (5  4) 

Thus  (5.4)  is  a  necessary  condition  for  the  variable-step,  variable-coefficient  .*I-Mabilit\ 
property  studied  here.  This  condition  is  satisfied  iff  (3  3))  holds.  Assume  that  tin*  coefficients 
are  smooth  enough  in  e  so  that  wc  can  write  fi(f)  =  —  rv(0)  4-  O(f').  Then  r )  >  0  is 
obviously  satisfied.  Furthermore,  for  arbitrary  f,  y(r)  >  0  iff  v(r)  >  t ,  p;irt  of  the  condition 
(4.1)  for  contracti vity  at  q  =  0.  Also,  for  |  e  |  <<  1,  <>(>  )  =  w(0)  —  rv(0)  +  and 

0(f)  >  —  is  equivalent  to  (4.7),  one  of  the  constraints  for  contractivitv  at  q  =  *.  Finally,  the 

4 

condition  y(r)  <  4<$(f)0(f)  becomes  -e  +  v(0)  <  v(0)-2f v2(0)-2fh*(0)  4-  <?U‘)  or ,  to 
0(f),  -e  <-2fv2(0)-2rw(0)  which  is  true  for  any  r,  I  r  |  <<  l,  iff 

H’(0)  =  - —  v2(0)  (5.5) 

2 

holds.  Condition  (5.5),  together  with  (4.7)  which  for  r  =  0  is  w  (0)  >  -  - ,  implies  v(0)  <  . 

4 

and  y(0)  >  0  requires  v(0)  >  0.  Thus  altogether,  for  |  r  |  <<  1  the  necessary  condition  (5  4) 
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for  ^-stability  is  satisfied  iff  both  (5  5)  and 


0  <  r<(>)  < 


(  5  M 


hold  which  for  f  =  0  are  equivalent  to  (3.15a)t  and  (3.15b).  respectively,  defining  the 
/4-contractive  formulas 

For  arbitrary  f  ^  0,  the  formal  algebraic  ^-stability  conditions  (3.‘>)  yield  the 
constraints  derived  hereafter.  The  conditions  <$  >  0,  .p  >14,  and  y  >  0  yield 


<  I  +  r*  >  n 

v  - .  r  0. 


(5  7) 


H-  >  +  ft*. 


(5.8) 


and 


V  >  f , 


(5  «) 


respectively.  For  any  f.  the  constraint  y  <  4<5tt  takes  the  form  r«I»(v.  k\  f)  >0.  where 


<l>(  V.  H’,  f)  :=  2|2fi-(l  +  c  *" )  |  h’  —  2  ( I  +  rV  +  r(3  +  r)v  + 


<5  10) 


It  is  satisfied  iff  ^(v,  w,  r)  -  0,  f>  0.  The  equality  4>(v,  w.  t )  =  0  can  be  written  in  the 

<  < 

form 


K-  =  [2(1  +  f 2 )  »'2  -  f  ( 3  +  e2)v-(l-r2)H2l2ff-(l  +  f2>1|  ~ 


(5.11) 


2  4  2  2  / 

a  hyperbola  whose  center  lies  at  h*  =  (2  +  e  +  e  )/4f  ,  v  =  (1  +  t  )/2 r  and  whose  two 
asymptotes  are  v  =  ( 1  +  f")/2e  and  w  *  [(1  +  f2)/2f\y  -f  [( 1  — r2)/4f " |.  This  hyperbola 
passes  through  the  points  ^  w  »  ^  ^  ,  v  =  ^  and  v  ~  rcPrcscnl“ 

ing  the  Trapezoidal  Rule  (77?)  and  the  two-step  Trapezoidal  Rule  (2 77?),  respectively. 

The  hyperbola  subdivides  the  (h\  v)-plane  into  a  connected  set  containing  the  minor 
axis  and  into  a  disconnected  set  consisting  of  two  branch  sets  and  containing  that  part  of  the 
major  axis  lying  beyond  the  vertices.  For  c  >  0,  the  "good"  region  defined  by  <J>  >  0  is  the 
disconnected  set.  Because  the  first  asymptote  is  the  borderline  case  of  condition  (5.7),  only 
one  of  the  branch  sets  is  "allowed"  by  condition  (5.7)  and  all  points  of  that  set  do  satisfy  that 
condition.  Consequently,  for  r  >  0  the  set  of  all  points  satisfying  the  formal  /f-stability 
conditions  is  the  (bounded)  intersection  of  (5.8),  (5.9)  and  >  0  It  is  represented  (for 

e  *  +  1/2)  by  the  shaded  region  of  Fig.  3  which  contains  the  arc  of  all  ^-contractive 
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methods.  For  comparison  with  the  other  figures,  the  i40-conlractivity  triangle  is  also  shown 
For  f  <  0.  the  "good"  region  <t>  <  0  is  the  connected  set  The  constraint  (5.7)  is  implied  by 
(5.9)  and  thus  the  active  constraints  are  again  the  same  For  r  =  -1/2,  their  (unbounded) 
intersection  is  represented  by  the  shaded  portion  of  Fig.  4  which  again  contains  the 
/4-contractivity  arc. 

Figures  3  and  4  indicate  that,  for  f  near  zero,  the  arc  of  the  A -contractive  formulas 
nearly  coincides  with  the  boundaries  of  the  shaded  regions  for  f  >  0.  resp  f  <  0.  Hence 
assuming  smoothness  for  f  near  zero,  only  the  methods  along  the  arc  are  admissible,  as  we 
concluded  above.  However,  wc  have  not  shown  that  a  non-smooth  strategy  would  not  work. 

Note  that  the  BDF  .  which  for  f  =  +  1/2  is  (9/S.  3/2)  anil  for  f  =  -1/2  is  (1/X. 
1/2),  lies  in  the  "stable  set"  only  for  f  <  0.  ft  was  mentioned  in  (19|  (hat  for  the  test 
problem  defined  by  (1.12)  and  (5.l)-(5.3)  with  X ( r )  <  0,  /testability  of  the  BDF  was  lost  for 
r  >  1  +  sf~2  (r  >  v^2/(2  +  y/~2)).  Indications  for  instability  at  smaller  step  ratios  (r  >  1.2) 
were  given  in  [1]  for  v  =  Ax.  A  constant,  complex,  in  the  sense  that  for  such  ratios  a  sufficient 
condition  for  contractivity  in  a  polygonal  norm  was  violated.  Here  it  is  show  n  that  for  some 
complex  A(/)  (oscillatory  problems)  the  BDF  actually  becomes  destabilized  by  the  problem 
defined  by  (1.12),  (5.1)  -  (5.3)  with  r  =  1  +  f  for  any  f  >  0. 


s 
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Figure  2.  ^-contractive  and  /4-contraeiivc  formulas  for  decreasing  steps 
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8  > 0 v <  5/4  BDF  (9/8, 3/2) 
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Figure  3.  ^-stable  ami  /I -contractive  formulas  for  increasing  steps 


Figure  4.  /(-stable  and  /(-contractive  formulas  for  decreasing  steps. 
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